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Abstract A recipe is presented for interpreting non-invasively the transport 
processes at work during relaxation of a cylindrical, supcrfiuid-fillcd vessel, 
after it is accelerated impulsively and then allowed to respond to the viscous 
torque exerted by the contained fluid. The recipe exploits a recently published 
analytic solution for Ekman pumping in a two-component superfluid, which 
treats the back-reaction self-consistently in arbitrary geometry for the first 
time. The applicability of the recipe to He II, 3 He, 3 He- 4 He mixtures and 
Bosc-Einstcin condensates is assessed, and the effects of turbulence discussed. 
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1 Introduction 

The study of transport processes in superfluids has a long history [TJ[3J[3]. 
Measurements of transport coefficients, like the viscosity or mutual friction, 
help to shed light on microscopic physics like the phonon- vortex and phonon- 
roton interactions in He II [1] . Transport processes are also useful for charac- 
terizing superfluids with complicated phase diagrams like 3 He [5] and 3 He- 
4 He mixtures [BJ, changing more strongly between phases than some thermo- 
dynamic variables like density. In non-laminar flows, the effective viscosity 
offers a practical way to quantify complex dynamics like vortex tangles and 
turbulence ^EM- 

Spin-up experiments have played a crucial role historically in the under- 
standing of superfluidity and represent one important arena where transport 
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processes play out [T0l[TTl[T2 II 1 3 IfM] . A diverse range of phenomena were ob- 
served and explored by Tsakadze and Tsakadze [T5], including spin down of 
smooth- and rough-walled containers, quasi-periodic oscillations, and spas- 
modic jumps in angular velocity. The data have never been reproduced quan- 
titatively by a first-principles theoretical analysis; in a recent publication [16] . 
we demonstrated agreement at the 0.5% level between the data and a hydro- 
dynamic theory for the smooth spin down, but the other phenomena remain 
unexplained 0- Subsequent experiments in helium II and 3 He have focused on 
quantum turbulence [18,19,20,21]. Below a critical temperature, the energy 
flux in the turbulent (Kolmogorov) cascade transfers to individual vortex 
lines in a Kelvin- wave cascade when it reaches the quantum scale, and the 
viscosity drops dramatically [9,22, 23 . Mutual friction and viscosity also play 
an important role in the onset of turbulence in 3 Hc- 4 Hc mixtures, which have 
been studied less thoroughly to date pMll2"5"] . 

In a spin- up experiment, the container responds to the viscous torque 
exerted by Ekman pumping |26j . The latter process is sensitive to transport 
physics that cannot be probed in traditional, small-amplitude, second-sound 
or torsional-oscillator measurements of the mutual friction and viscosity. For 
example, viscosity measurements using fine-wire vibration are difficult to 
conduct in 3 He, where the oscillations are strongly damped. Ekman pumping 
can also be used to understand turbulent flow, where the spin-up times are 
typically faster than laminar flow [27,22,28]. In astrophysical applications 
like neutron stars, remote observations are the only possible way to measure 
the supcrfluid coefficients in bulk nuclear matter |29U30] . Spin- up studies are 
independent of experiments involving vibrational modes. 

The paper is structured as follows. In iJ5J we outline the set-up of a generic 
spin-up experiment and the role potentially played by turbulence therein. In 
321 an analytic solution for supcrfluid Ekman pumping derived by van Eysdcn 
and Melatos [3"TU3"2"] is used to demonstrate how to interpret the measured 
container response in terms of elements of the flow physics and extract ac- 
curate values for the transport coefficients. In we discuss how Ekman 
pumping experiments can be conducted and interpreted in He II, 3 He, 3 He- 
4 Hc mixtures and Bosc-Einstcin condensates, with an emphasis on identifying 
the optimal regimes in which to probe various transport phenomena. 



2 Spin-up experiment 

Consider a rigid, cylindrical vessel filled with supcrfluid. To begin with, the 
vessel is spun up to a constant angular velocity ft s about the z-axis (unit 
vector k), until corotation between all components is established. The angular 
velocity of the vessel is then increased impulsively by a small amount 8fl. 
The subsequent motion of the vessel is tracked as accurately as possible, as 
it responds to the viscous torque exerted by the superfluid and the friction in 
the mount. The angular velocity versus time is fitted to the spin-down curve 
predicted by theory [3~1.32 ■ Two exponential time-scales arc identified, from 



1 Adams et al. [T7] explained the spin up from rest in rough-walled containers 
below 1.3 K by modelling pinning as friction between the vortices and the walls. 
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which the Hall-Vinen mutual friction parameters B, B' and shear viscosity 
r} can be extracted. 

To analyze and interpret the post-jump response, we exploit a recently de- 
rived theoretical solution for the spin up of a superfluid-fillcd container [3TJ 
[32] based on the Hall-Vinen-Bekharevich-Khalatnikov (HVBK) equations, 
which model the superfluid in terms of viscous and inviscid components cou- 
pled by a mutual friction force. The analysis assumes that the flow is linear, 
with Rossby number e — 5Q/Q S <C 1. Linearized and in the rotating frame, 
the HVBK equations take the dimensionless form [331131] 



E l/2^n +2kxVn = _ Vpn + EV 2 Vn+ P lF (1) 
OT p 

£i/2^i + 2 kxv s = -V Ps -^F, (2) 

OT p 

V • v n = , (3) 

V • v s = . (4) 

where v„ jS and p n . s are the bulk velocities and pressures of the viscous (n) 
and inviscid (s) components. In the theory, it is assumed that the container 
rotates with Ekman number £<Cl. The Ekman number is defined as 

(5) 



h 2 f2. 



where v n = r\j ' p n is the kinematic viscosity of the viscous component. The 
linearized mutual friction force takes the Hall-Vinen form [34] 

F = Bk x [k x (v„ - v a )] + B'k x (v n - v.) , (6) 

and B, B' symbolize the mutual friction coefficients. The velocity and pres- 
sure scales are chosen to be h8f2 s and pfi s h 2 Sf2 s respectively, where p is the 
total mass density of the fluid, and h is the half-height of the cylinder. 

In {l}-©, both fluid components are assumed to be incompressible and 
the thermal conductivity vanishes. Therefore the only dissipative coefficient 
appearing in H])-© is the first viscosity v n [35]. The circulation in the invis- 
cid component arises from a rectilinear array of quantized vortices polarized 
along the rotation axis, whose critical angular velocity for formation is [36] 



where r is the vessel radius, ao is the vortex core radius, and r = 2i:h/m 
is the quantum of circulation (Planck's constant divided by the mass of one 
superfluid boson). The HVBK equations are derived by averaging smoothly 
over a high density of vortex lines and require the vortex line separation to 
satisfy (P/2i7 s ) 1 ^ 2 <C r [3"71I3"51I5M3"5] . The inviscid component spins up when 
vortices created at the side- walls of the cylinder migrate inwards [ID] . The 
new vortices added during the impulsive spin up are a small fraction of the 
total number (oc e), so the analysis is valid even in the presence of a vortex 
formation barrier 1211. In other words, even if new vortices cannot be added 
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during spin-up and thus a vortex-free region forms close to the cylinder wall, 
the overall dynamics would not be significantly affected, since the vortex-free 
region is small. 

To keep the problem analytically tractable, the walls of the container are 
assumed to be smooth and pinning is neglected, as contrived in many labo- 
ratory experiments, e.g., Refs. j"PWTT" in . The vortex tension, parameterized 
by v s , is also neglected □■ We do this, not because the tension is necessarily 
small, but rather because v s only appears in the inviscid boundary layer solu- 
tion [39] ; the interior flow solution and spin- up time of Ref. [39] are recovered 
despite assuming v s = 0, suggesting that v s plays a minimal role in the spin- 
up of smooth-walled containers. (See section 4.4 of Ref. [3T] and section 2.4 
in Ref. [TB] for details.) Indeed, the v s = theory agrees well with experi- 
mental data for He II in smooth-walled containers in the temperature range 
I.4K<T< 1.8 K [T5]. In the present analysis, the viscous component obeys 
the usual no-slip and no-penetration boundary conditions with respect to the 
tangential and normal velocity respectively For the inviscid component, no 
penetration is sufficient. The theory is not valid in rough-walled containers, 
where the spin up of superfluids proceeds differently [THTllll"Tll3"9"] . 

Angular momentum is transferred from the container to the fluid by Ek- 
man pumping, which is illustrated in Figurc[T] The flow pattern of the viscous 
(left panel) and inviscid (right panel) components is shown at t = E^ 1 / 2 flj 1 . 
In the viscous component, the Coriolis force drives a secondary flow in bound- 




Fig. 1 Ekman pumping for a two-component superfluid in cylindrical geometry 
at t = E~ 1 ^ a f2^' 1 for the viscous component (left panel) and the inviscid compo- 
nent (right panel). Poloidal streamlines are plotted as thick black curves. The flow 
is clockwise in the top-left corner and mirrored about the vertical and horizontal 
bisectors in the other quadrants. The magnitude of the azimuthal velocity is rep- 
resented by the shading, where red is highest and blue is lowest. Model parameters 
are B = 0.8, B' = 0.2, p n = 0.5/0, E = 0.01, K = 0.25 and T = -0.0063. 



2 The tension parameter v s looks like a kinematic viscosity, but it is non- 
dissipative and gives rise to Kelvin waves. It takes the form v s = (r/An) ln(6o/ao), 
where bo = (F/2f2 a ) 1 ^ 2 is the inter- vortex spacing. 
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ary layers of thickness E 1 l 2 h along the top and bottom of the cylinder. The 
boundary flow draws in fluid from the interior, spins it up, and cycles it back 
into the interior on the Ekman time-scale E" 1 ! 2 fl~ x . The inviscid component 
is spun up by the viscous component via the mutual friction force. When the 
mutual friction is strong, i.e. B 3> E 1 / 2 , as in Figure [TJ the inviscid com- 
ponent undergoes poloidal circulation as it is dragged along by the viscous 
component. For intermediate mutual friction coupling, i.e., B ~ E 1 / 2 , the 
inviscid component does not circulate poloidally, and the azimuthal viscous 
component spins up its inviscid counterpart over the mutual friction time- 
scale. For weak coupling, i.e., B -C E 1 / 2 , the viscous component spins up 
uninhibited by the inviscid component, which is spun up by mutual friction 
over a much longer time scale. A detailed study of the role of mutual friction 
in the spin- up of two-component supcrfluids can be found elsewhere [31] . 

The Ekman and Rossby numbers control a number of important aspects 
of the experiment. Both should be made as small as is practical, although for 
E < e 2 turbulence may switch on, with important consequences (see below). 
The scale of the apparatus h enters through the Ekman number; a larger con- 
tainer reduces E. For a container with h = 2 cm rotating at fi s = lrads -1 
(as in Ref. [IS]); we obtain E <C 1 for most known terrestrial supcrfluids 
(or indeed fluids), including helium II. Before the impulsive spin up, a con- 
stant rotation rate should be maintained for a few minutes (corresponding 
to several Ekman times; the diffusion time is hours), to ensure that the in- 
terior fluid attains steady-state co-rotation and all inertial oscillations are 
damped. The Ekman number also controls the maximum duration of the 
spin-up event, which one contrives to be shorter than the Ekman time-scale 
E~ X I 2 Q~ X [16]. The Rossby number must be large enough for efl s to be 
accurately measurable. 



2.1 Turbulence 

The above analysis also describes turbulent flow, when the transport coeffi- 
cients are interpreted as eddy-averaged quantities. Many qualitative aspects 
of laminar Ekman pumping persist when the flow is turbulent; e.g., the sec- 
ondary flow pictured in Fig. [1] is still present [41]. Similar qualitative be- 
haviour has also been observed for the transition to turbulent vortex flow in 

He II pang. 

In a superfluid, turbulence is either quasi-classical, quantum mechanical, 
or a combination of the two. Quasi-classical turbulence occurs on macroscopic 
scales through the formation of Kolmogorov-like eddies, even when there 
is no counterflow (i.e., v„ = v s everywhere). Quasi-classical turbulence is 
predicted to occur when 

Re a = ^^, (8) 

which is the ratio of the reactive and dissipative components of the mutual 
friction force, becomes large [421127] , Equation (JSJ suggests that this occurs 
when p n — > as T — > 0. Turbulence can also enter classical Ekman pumping, 



6 



when the Ekman layer becomes unstable [4111391116] . In a Navier-Stokes fluid, 
this occurs for 

£ < e 2 (56.3 + 58.4er 2 , (9) 
suggesting that e < E 1 / 2 <C 1 is required for laminar Ekman pumping. 

Quantum turbulence arises when the counter-flow velocity along the vor- 
tex lines, (v„ — v s ) • fi, exceeds the critical value \vdg\ = 
at which point Kelvin waves are excited and a fully developed, self-sustaining, 
reconnecting vortex tangle develops |^1[2U1I^T1H5] . This type of turbulence 
occurs on microscopic scales (that of the vortex lines), even if v„ and v s 
are uniform on macroscopic scales. The mutual friction force per unit mass 
(written in terms of starred, dimensional quantities) takes the Gortcr-Mcllink 
form USES] 

F* = A ^(v:-<) 3 , (io) 

where A is a phcnomenological constant which depends on the type of fluid 
and temperature. In Ekman pumping, the counter-flow velocity scales as 
eE 1 / 2 Q s h [21], so that we require 

e<2{2 Vs /v n ) 1/2 (11) 

for laminar flow. We check this condition in several applications in §Q 

To characterize the type of turbulence, the superfluid Reynolds number 
is defined [37] 

OR 2 

Re s = — . (12) 

Turbulence is observed in He II and 3 He-B when both Re Q > 1 and Re s 3> 
1 hold [UJ. Quasi-classical Kolmogorov turbulence is found in the regime 
Re a < Re*/ 2 , while quantum turbulence is found when Re a > Re 1 / 2 . The 
turbulent state depends on the particular fluid in question; we discuss some 
particular examples in 21 

An important issue is how the dissipation changes between laminar and 
turbulent flow and how this carries through to the observed spin-down time of 
the container. Typically, turbulence enhances disspiation. This is observed in 
spin-down and spin- up experiments in 3 He-B, where the turbulent response 
is faster at a given temperature [22 , 27. 28 . Interestingly, experiments in He II 
show that the spin-down time is longer than expected for laminar flow near 
the superfluid transition temperature, where turbulence is predicted |16j . 
Dissipation can also change significantly within the turbulent state. Below 
a critical temperature, the energy flux in a turbulent (Kolmogorov) cascade 
transfers to individual vortex lines in a Kelvin-wave cascade when it reaches 
the quantum scale, and the viscosity drops dramatically [9"l [2"2"ll2"5] . 



3 Container response 

3.1 Spin-down curve 

Let fi{r) = fi s + 5f2f(r) denote the angular velocity of the vessel at r > 0, 
where the dimensionlcss variable r is the usual time-coordinate normalized 
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by the Ekman time E^^Qj 1 [261131"] . At time r = 0, when the cylinder 
is released after being accelerated impulsively, the initial velocities of the 
viscous and inviscid components are zero in the rotating frame, and f(0) 
equals unity. The constant frictional torque in the apparatus is denoted To 
in dimcnsionless units and is negative. For a rigid, symmetric container, the 
integral equation of motion is |32| 



/ (r) = -K f dr' [g A (r - r') + g B (r - r')} f(r') + T r + 1 
Jo 



(13) 



Equation (jT3"]) captures self-consistcntly the back-reaction of the supcrfluid 
on the container and vice versa, e.g. by including the continuously changing 
boundary conditions on the superfluid, as Q{t) evolves. For a cylinder of 
radius r, height 2h, and moment of inertia I c , the ratio of the nominal moment 
of inertia of the fluid (rotating as if it were a rigid body of uniform density 
p) to the moment of inertia of the container is given by 



K 



pwr^h 



(14) 



In JT 



we make the definitions 

J (e"+ T 



9 A (r) 
9 B (r) 



LO + — OJ- 

LJ _ ( e "+r _ i) _ w+ ( e u-r _ 



(15) 
(16) 



uj + — UJ- 

where the overdot symbolizes a derivative with respect to r, and we have 

1/2 



U± 



I 



± 



f3J 



2BE- 1 ' 2 
2-B' 
2X+ 
X\ + A 2 _ 5 
J = p n \- , 



(17) 

(18) 

(19) 
(20) 



A+ = 



(B' -2) 2 + B 2 
{p n B' - 2f + (p n B) 2 



1/2 



-F 



2 Ps B 



(p n B> - 2) 2 + (p n B) 2 



Equation (fT5)l has the exact solution 

/(r) = (1 - C - exp {w+t) + Ccxp (w.t) + Dt + f 



1/2 



(21) 



(22) 
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Fig. 2 Spin down of a superfluid-filled cont ainer accelerated impulsively at r = 0. 
The thick black curve corresponds to ()22p . while the first and second terms in 
(12211 are graphed separately as blue-dashed and red-dotted curves respectively. The 
model parameters are the same as in Figure [T] 



with the auxiliary definitions 
'0 + 1 + JK 

w± 



± 



+ 1 + JK 



0J (1 + K) 



1/2 



c 

D 



Kw + (m- + 0) (w- + T ) 
(1 + K) -w-) 
To 



1 + K 
1 

1 + K 



1 + To 



VJ- + ZU+ LJ- + UJ + 



, (23) 
(24) 
(25) 
(26) 



Equation (|2"2"j) can be verified by substitution into ([TBI . It comprises two 
exponential decays, which depend on the mutual friction, viscosity, and the 
relative inertia of the container; a linear spin down, in response to the con- 
stant external torque; and a steady-state solution, which the system would 
relax to in the absence of an external torque. The solution (f2"2"j) is graphed as a 
thick black curve in Figure [5J The two exponential decaying terms in (|22[) are 
also graphed separately as blue-dashed and red-dotted curves respectively. 
For these parameters, both decays are clearly and separately evident in the 
overall recovery 

The detailed form of the spin down depends on the strength of the mu- 
tual friction; see Ref. [31] for a detailed discussion. When the superfluidity 
coefficients satisfy |3 > 1, the two fluid components are drawn together by 
mutual friction much faster than the Ekman time. In this limit, we obtain 
-n- 1 E- 1 l' 2 w^ 1 = E^^n^Jil + K)/h, which is the superfluid Ekman 
time [3"§ll3"T] . while — Q~ 1 E^ 1 / 2 wZ 1 = Bfl^ 1 is the coupling time-scale for 
the two components due to mutual friction. We also obtain C 1 — C — /oo, 
and the spin down is dominated by the first term in (|22"]) . This strong cou- 
pling limit is discussed further in M3.31 For intermediate (i.e., ~ 1) coupling, 
both exponentials in (f2"2"j) are important, and — Q' 1 E~ 1 ^ 2 zu^. 1 are complex 
mixtures of the Ekman and mutual friction coupling times. 

The spin-down time of the container is strongly dependent on K. Increas- 
ing I Cl e.g. by attaching a metal disk to the spindle on which the container 
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rotates, reduces K and increases the spin-down time. This may be desirable 
if the container response is too rapid for the shape of the spin-down curve 
to be measured accurately. As an independent check, the experiment can be 
repeated for different values of K. 

The solution (|22|) assumes that the viscous and inviscid fluid components 
are initially co-rotating. As a consequence, the pre-factors 1 — C — foo and C 
are always positive, and the spin down is monotonic. In general, either 1 — C — 
/oo or C can be negative, when the inviscid and viscous components initially 
rotate differentially. This can cause an overshoot: the angular velocity of the 
container drops below the steady-state value before rising again. Although 
this possibility is interesting, especially in astrophysical applications [3D], we 
do not consider it here, because the initial velocities of the HVBK components 
are difficult to fine-tune experimentally. 

Dual exponential recovery is characteristic of the above system, irrespec- 
tive of geometry. Qualitatively similar behaviour is seen in a sphere, ex- 
cept that the shape of the spin-down curve is quasi-exponential (strictly, 
a latitude- weighted sum of exponentials) |31j . The hydrodynamic torque in 
non-parallel-plate geometries varies along the boundary, and solutions to (fTB"]) 
cannot be written in analytic form like (|22p . The solution governing equations 
for /(t) in general geometry is given in the Appendix. 



3.2 Extracting B, B' and E from the spin-down curve 



The analytic solution (|22p can be inverted, so that the supcrfluid transport 
coefficients B, B' and E can be inferred from the experimentally measured 
spin-down curve Q(t). From (|22"]) - (|26|) . we obtain 



T = 

W+VJ— 
ZU+ + n7_ 



—KJ + T , 

0J(l + K), 
-((3 + 1 + JK) 



(27) 

(28) 
(29) 



Restoring dimensions according to uj± = zu±Q s E 1 ' 2 , t = E 1 ' 2 O s V, and 
T * = Snn s E^ 2 I c T , equations ([27 | - (|29)l can be rearranged to give 



JE 1/2 = [5f2{2 s K 
/3E 1/2 = w\vj*_ 
IE 112 = -vo*vo 



1 


1 


fdn\ 




Jc ' 





K5f2 



[l + K)Q s 
K8Q 



T * 



Ic 



(1 + K) Q s 



Tl 

Ic 



dn 
It 

~dt 



t=0. 



dsn 

~dt 



a 



(30) 
(31) 

(32) 



The left-hand sides of equations ([30|) - ((3"2"|) feature the unknowns B, B' and 
E, while the right-hand sides feature measurable quantities. For example, 
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when the container is empty, its spin-down rate is dfl/dt = Tq /I c , giving the 
external torque in the apparatus. When the container is full, we require the 
container to accelerate slowly, with Tq/I c (1 + K)Q 2 S E X I 2 C 1 to satisfy the 
linearity assumption. 




012345 12345 



B B 

Fig. 3 Contours of /3E 1/2 (dashed curves) and I/J (solid curves) for p„ = 0.1 (left 
panel) and p n = 0.9 (right panel). The mutual friction coefficients B and B' can 

be read off easily from the intersection of the fiE 1 ^ 2 and I/J contours, measured 
from the spin-down curve. 



The quantities (3E 1 / 2 and I/J depend only on B, B' and p n . If p n is 
known by some other means (e.g. published data from previous experiments, 
such as Ref. [IE]), equations (j3Tj) and (|32"j) divided by ([3D]) are sufficient to 
extract B and B' . This is done graphically in Figure [3] Contours of constant 
I/J arc plotted as solid black curves on the B-B' plane, while ftE 1 / 2 contours 
are plotted as dashed curves on the same plane. Values of // J and /3E 1 / 2 are 
calculated from equations (I3D1) (f35|) . Given measurements of the quantities 
in the right-hand sides of ([50]) (|32|), B and B' can be read off the vertical 
and horizontal axes at the intersection point of the curves. Once B and B' 
have been found, E is determined uniquely from (|30[) or (|32j) . 

Cylindrical geometry is used here as an example, because a concise ana- 
lytic solution can be obtained. Solutions in general geometry also exist, but 
they do not simplify like in (|3D]l or ([32"]) , so the fit must be done by adjusting 
7], B and B 1 by trial and error [3D]. The solution to (f2"2"j) in general geometry 
is presented in the Appendix. 



3.3 Strong coupling 

Many supcrfluids are strongly coupled by mutual friction, in the sense that 
one has /3 > 1. For example, He II has 1.526 < B < 12.693 for 1.30 K < T < 
2.1767 K giving f3 ^> 1 in the Ekman pumping regime £ « 1 [49]. In this 
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limit, (fT3j) reduces to 



f(r) = (l+Ky 



+T t + 1 + 



K 



1 - 



J (l + ir) 



exp (n7 + r) 



TnK 



J (1 + K) 



(33) 



With one of the exponential decays now imperceptible, there are fewer mea- 
surable quantities available to constrain B, B' and E [unless the cxp(n7_r) 
term can be resolved by a sufficiently sensitive experiment]. Equation (I33P 
does not contain I and (3; all the dependence on B and B' resides in J. 
Restoring dimensional variables, there is now only one relation between B, 
B' and E, namely (J3DJ). 



4 Applications 

4.1 Helium II 

Helium II is the most extensively studied superfhhd. Its transport coeffi- 
cients have been measured and independently verified in several ways. Table 
[T] summarizes a representative sample, noting the experimental technique 
and temperature range T. The entries draw on information published elec- 
tronically by Donnelly |49| . Historically, viscosity measurements have been 
attempted with rotating cylinders [501151] . damped fine- wire vibration |521 
[S3], torsional quartz crystals [SI], and torsional oscillators [SSJ. The rotat- 
ing cylinder viscometer measures the shear viscosity rj directly, whereas the 
other methods measure the kinematic viscosity r)/p n , so that p n must be in- 
dependently measured to extract rj. The methods agree within experimental 
errors. With regard to mutual friction, the B coefficient is determined from 
the excess attenuation of second sound due to the presence of vortices in a 
rotating container 56,57,58,59,60 , while B' is measured by the removal of 
a resonant degeneracy in a rotating second-sound cavity [61,58,59 ; see Rcf. 
[r?2"] for a review. Values for B and B' were compiled by Barenghi et al. [S2] 
and are also available electronically in Ref. |49j . 

The Ekman technique, described in SI3 - 2 1 for extracting rj, B, B' from spin- 
down curves cannot compete with the techniques in Tabled] for helium II, as 
it is less accurate, and spin-up experiments are more finicky. However it does 
possess two virtues. First, it is physically independent of the other techniques, 
relying on Ekman pumping instead of vibrational modes, and can therefore 
serve as a cross-check in new situations where the data are still uncertain. 
Second, it is sensitive to nonlinear phenomena like turbulence, which has been 
studied recently in experiments [H, 19, 20,21] and is usually (intentionally) 
not excited in experiments involving small-amplitude vibrational modes. We 
enlarge on this second point below. 

At low temperatures, where the viscous fraction is small, one has Re a ~ 
10 3 and turbulence is ubiquitous in He II [5]. Below a critical temperature, 
the energy flux in Kolmogorov cascade reaches down to the quantum scale 
and a Kelvin-wave cascade forms, dropping the effective viscosity by an order 
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Table 1 Measurements of transport coefficients in helium II 



Method 


Quantity 


T(K) 


Refs. 


Oscillating disks 


'/ 


(l.Z, I.U) 


uasn lay lor |4o] 


Couette viscometer 


'1 


(0.78, 2.079) 


Woods & Hollis Hallett [50] 


Vibrating wire 


V/Pn 


(1.52,2.16) 


Tough et al. [52] 


Vibrating wire 


V/Pn 


(1.2,4.2) 


Goodwin [53] 


Quartz crystal 


V/Pn 


(1.75,2.195) 


Webeler & Allen [5J 


Torsional oscillators 


V/Pn 


(1.8,4.4) 


Wang et al. [2>] 


Quartz tuning fork 


v/pn 


(0.1,2.2) 


Zadorozhko [I] 


Second sound 


B 


(1.289,2.060) 


Hall & Vinen [55] 


Second sound 


B' 


(1.28,2.156) 


Snyder & Linekin [61] 


Second sound 


B 


(1.2,2.08) 


Bendt [57] 


Second sound 


B, B' 


(1.383,2.167) 


Lucas [5S] 


Second sound 


B, B' 


(1.75,2.167) 


Mathieu et al. [59] 


Second sound 


B 


(1.40,2.05) 


Miller et al. [BD]^ 



of magnitude [9]. For T < 0.5 K, the viscous component also becomes a gas 
of ballistic quasi-particles, and the HVBK description no longer applies [B2]- 

As the temperature approaches the lambda point, v n decreases. Depend- 
ing on e, fl s and h, the condition Q may be met. For the experiments of 
Tsakadze and Tsakadze, where h = 2 cm, Q s ~ 3rads _1 and e ~ 0.15, 
turbulence is expected for T > 1.6 K [16]. Decreasing /i, l? s and e excites 
turbulence at higher temperatures. In He II, we have 0.5 < v s jv n < 6.1 for 
1.3K < T < 2. IK, varying approximately linearly in between [32]. Hence, 
the Glaberson-Donnclly instability (fTT| is not excited for the experiments in 
Ref. [15) . and a window exists where laminar flow is possible in He II. 

Recently, excellent agreement at the 0.5% level has been obtained be- 
tween the theory of Refs. [3T1I32] and the experiments of Tsakadze [15] in the 
temperature range 1.4K < T < 1.8K [TB]. Agreement is obtained using the 
generally accepted values of r], B and B' from Ref. [13] ■ This is strong evi- 
dence that the flow is laminar in this temperature range, and that the HVBK 
equations accurately describe the flow. In this way, the spin-down analysis in 
32] complements the traditional approaches in Table Q] by indirectly probing 
the conditions for turbulence in such systems. 

In He II, the mutual friction is strong [~ 0(1)] [1PJ- Therefore, only 
one time-scale is observed, and B, B' and E cannot be inferred uniquely. 
However, if the kinematic viscosity is known from other experiments [49], 
we can place useful constraints on B and B' . We illustrate how in Figure 
3] with an arbitrary example. Suppose that, from a spin-down experiment 
in He II at T = 2.0 K, we measure JE 1 / 2 = 0.002 ± 0.0002 using ([3D]). We 
find E = 9.11 x 10~ 6 using Ref. [49] (corresponding to a container with 
h = 2 cm rotating with fi s = 5rads _1 ). The thick black curve corresponds 
to JE 1 /' 2 = 0.002, while the shaded region bounded by the dashed curves 
represents the range JE 1 ' 2 = 0.002 ± 0.0002 for E = 9.11 x 10~ 6 . The 
generally accepted values of B = 1.008 and B' = 0.043 at T = 2.0 K are 
marked as a black dot on the figure. 
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1 2 3 4 5 



S 

Fig. 4 Constraint on the mutual friction coefficients from a hypothetical measure- 
ment. Contours of JE 1/2 = 0.002 ± 0.0002 for E = 9.11 x 10" 6 (shaded region 
between curves) and J_E 1//2 = 0.00 2 (t hick black curve). The values of B and B' 
from the adopted database in Ref. [49] are plotted as a black dot. 



Table 2 Transport coefficient measurements in 3 He 



Method 


Phase 


Quantity 


T (mK) 


Bo (T) 


Refs. 


Vibrating wire 


A, B 


V/Pn 


(1,100) 


(0,0.9) 


Alvesalo et al. [||] 


Torsional oscillator 


A, B 


r)/Pn 


(1,100) 


(0,0.9) 


Alvesalo et al. [54] 


Vibrating wire 


B 


V/Pn 


(0.6,2.49) 





Carless et al. [55] 


Torsional oscillator 


A. B 


V/Pn 


(2.11,2.64) 





Main et al. [66] 


First sound 


B 


V 


(1.2,2.64) 





Eska et al. [67] 


Torsional oscillator 


A 


nip™ 


(2.1,2.64) 





Hook EE] 


Vibrating wire 


A 


ri/pn 


(1.85,2.64) 


(0.5,9.2) 


Hata et al. [69] 


Torsional oscillator 


A 


V/Pn 


(0.3,3) 


(0.6, 12) 


Akimoto et al. |70| 


Torsional oscillator 


B 


V/Pn 


(0,5) 





Nakagawa et al [71| 


Torsional oscillator 


A 


V/Pn 


(1,80) 


(0.5,14.6) 


Roobol et al. |72j 


Kapton diaphragm 


A, B 


B, B' 


(0.5,2.49) 





Bevan et al. |73] 


Time of flight 


B 


B 


(1,2.49) 





Finne et al. [71] 



4.2 3 He-B 

Measuring transport coefficients in supcrfluid 3 He is additionally problem- 
atic because of the complicated texture physics, arising from spin interactions 
between 3 He atoms [75] • Supcrfluid 3 He can exist in an anisotropic A phase 
or an isotropic B phase. Transport coefficients in the anisotropic A phase 
are usually formulated as tensor quantities; different tensor components are 
measured in different experimental situations. Table [3J summarizes a repre- 
sentative sample of viscosity and mutual friction measurements, noting the 
material, phase, experimental technique, temperature range T, and magnetic 
field strength Bq. The different methods agree within experimental errors 
[72]. 

In 3 He, the mutual friction cannot be measured using second sound, be- 
cause acoustic waves are overdamped by the high viscosity [73ll74| . Alterna- 
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tive approaches include measuring the normal modes of transverse vibration 
of a Kapton diaphragm, which is used to excite a superfiow |76j . and mea- 
suring the time of flight of a vortex injected into a rotating cylindrical vessel 
[74] . The resonant response of quartz tuning forks in 3 He has been investi- 
gated. As yet, however, tuning-fork experiments have not yielded accurate 
values for the viscosity or mutual friction |77IIT5] . 

In 3 He-B, turbulence is inhibited by a large viscosity pl2"ll21] . Noninvasive 
nuclear magnetic resonance measurements show that, in a cylinder, spin down 
is laminar all the way down to temperatures below 0.2 T c , despite Re a ~ 10 3 
[2"7] . For 3 Hc-B, one has v s ~ 10 -7 m 2 s -1 . Hence, from (fTTj) . we require 
e < 0.02 to prevent the onset of the Donnelly-Glaberson instability. Also, 
in 3 He, vortex cores are at least 10 2 times larger than in He II, and surface 
pinning and friction do not have a substantial impact on flows in smooth- 
walled containers j3^. However, like He II, the viscous component of 3 He 
becomes ballistic in the low-temperature limit and the HVBK description 
breaks down. 

For 3 He-B we have approximately 5 x 10~ 4 m 2 s~ 4 < v n < 2.5 x 10~ 4 m 2 s~ 4 
for 0.4 < T/T c < 0.8 [71], approximately 10 4 times larger than in He II. 
Therefore, we require fi s h 2 > 5 x 10~ 4 m 2 s~ 4 for Ekman pumping to occur 
(E <C 1). This condition is met in typical 3 Hc-B experiments. For example 
in Rcf. [401127] . one has h « 10cm, Q a w lrads" 1 , and hence E ~ 0.02. 
Therefore 3 He-B is amenable to the interpretative analysis described in SQ] 
and 33 

4.3 3 He- 4 He mixtures 

Equations ([I])-© describe general two-component fluid mixtures such as 3 He 
and 4 He just as well as an HVBK superfluid [751I501I3"?] . To date, the study 
of 3 He and 4 He mixtures has focused on their thermodynamic properties, 
for applications such as cryogenic refrigeration [5Tll82j . and the transition 
to turbulence [25] • However, transport coefficients and turbulence in 3 He- 
4 He mixtures have been studied to some extent in the temperature range 
0.15 K < T < 2.17 K, where dilute concentrations of 3 He in the normal liquid 
phase are added to superfluid 4 He. Castclijns et al. [H] investigated a 6.6% 
mixture of 3 Hc in 4 Hc for T < 150 mK and concluded that 3 He generates a 
vortex tangle in 4 He in the same way that the normal component generates a 
vortex tangle in the superfluid component in He II. The resulting mutual fric- 
tion has a Gorter-Mellink form (fTU)) . with A ~ 1.5 x 10 4 kg -1 ms, ^30 times 
larger than a vortex tangle in pure 4 Hc. In mixtures, 3 He also concentrates 
in the vortex cores of the 4 He, which may alter the dynamics significantly. 
Qin et al. [53] measured the shear viscosity of mixtures ranging from 1-5% 
3 He and observed that the shear viscosity decreases with increasing He con- 
centration. However, comparison with independent viscosity measurements 
using capillary flow j5HI55] and torsional oscillators [561151] yielded discrep- 
ancies of nearly 100% for 1% concentrations. Spin-up experiments sensitive 
to the independent Ekman process may help to resolve this controversy. 

Recent studies suggest that 5% concentrations of 3 He in 4 He enhance 
laminar flow stability and increase the critical velocity for the transition to 
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turbulence below 1 K [25] . This suggests that the spin-down experiments pro- 
posed in this paper apply to 3 He- 4 He mixtures. For example, the distinctive 
rise in the spin-down time of the container with the onset of turbulence could 
serve as a useful diagnostic. 

4.4 Bose-Einstein condensates 

Dilute-gas Bose-Einstein condensates (BECs) are an interesting application, 
as they are less noisily coupled to their environment than liquid helium [87] . 
Condensates rotating in magneto-optical traps form vortices which do not pin 
except when forced, e.g. by an optical lattice [HE]; traps therefore constitute 
smooth-walled containers. Another feature of the hydrodynamic analysis is 
that it can be applied to containers with arbitrary axisymmetric geometry, 
such as the parabolic traps used to construct BECs. The general theory is 
summarized in the Appendix. However, because the trap interacts with the 
whole BEC simultaneously and not just with its edges, no Ekman boundary 
layer forms. The theory in Ref. |31| still applies; one may imagine a situation 
in a two-component BEC where one component is spun up, leaving the other 
to respond via mutual friction. In a BEC, the vortex core size approaches the 
inter- vortex spacing [3T]. It is hard to create a rapidly rotating BEC with 
many vortices, because the condensate becomes unstable as the rotation fre- 
quency fi s approaches the trap oscillation frequency uj±. Also, as Q s — > 
the vortex core size grows until it saturates at the inter-vortex spacing and 
the vortices overlap [55]. Nevertheless, states with fl s w 0.995wx containing 
hundreds of vortices have been constructed successfully [ST), 91 , 92 . 93] , for 
which the HVBK approximation may be suitable [9T1I92] . 

4.5 Neutron supcrfluid: weak coupling 

Recently, much attention has been devoted to the thermodynamic phases and 
transport coefficients of bulk nuclear matter, which is generally supcrfluid 
and even superconducting [M] . This flowering of interest has its origin in the 
desire to test the theory of quantum chromodynamics in the MeV-energy, 
many-body regime. The hydrodynamics of nuclear matter in beta equilib- 
rium involves 17 transport coefficients in general [73] ■ The viscosity has been 
measured in the high-energy (GeV), few body (~400 nuclcons) regime by 
Au + Au collision experiments in relativistic heavy-ion colliders [951I5S] . The 
other coefficients have not been measured. Indeed, at present, the only prac- 
tical way to measure B and B' is by studying the phenomenon of rotational 
glitches in neutron stars [9"7l[3"0] . where Nature arranges a facsimile of the 
spin-up experiment discussed in this paper. Neutron star measurements of 
nuclear transport coefficients are used to investigate the existence of exotic 
nuclear supcrfluid phases like color-flavor locked and two-flavor superconduc- 

tors pii^i nroinui] . 

Weakly coupled two-component superfluids are believed to exist in the 
interiors of neutron stars, where an ultra-dense superconducting proton- 
electron plasma and a neutron condensate play the role of the viscous and 
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inviscid components respectively [102, 103, 104, 105,30 . Estimated densities 
range from 10 11 gem -3 near the surface to 10 14 gcm -3 or higher in the core. 
Theoretically, it is thought that the viscosity in the star arises from either 
neutron-neutron or electron-electron scattering, which give E k, 10~ 7 [1061 
1107] . while the mutual friction force arises from inelastic electron scattering 
off neutron vortex cores, typically giving 10~ 5 < B < 10~ 4 and B' ps B 2 [21)1 

The recovery of a neutron star after an impulsive spin-up event (or glitch) 
is believed to result from the hydrodynamic relaxation of the crust and fluid 
according to the Ekman pumping process analyzed in this paper |30j . In 
Figure [5j we plot the observed spin-down curve of the Vela neutron star after 
its 1985 glitch. The glitch recovery has been fitted to a function of the form 

una 

Q{t) =cie ro +* + C2e ro * t +c 3 < + C4, (34) 

with c 3 = to give a = (1.73 ± 0.01) x 10~ 5 , c 2 = (4.1 ± 0.6) x 10~ 7 , c 4 = 
(9.475 ± 0.006) x 10" 5 , w\ = (2.19 ± 0.05) x 10~ 7 , w*_ = (1.1 ± 0.1) x 
10- 5 , Q s = 70.4 (all in units rads" 1 ), and (df2/dt) t=Q = -8.69xl0" 7 rads- 2 . 
In a neutron star, K and p n are not known. The only restriction is that (|3"0"j) - 




50 100 150 200 250 



t (days) 

Fig. 5 Observed spin-down curve of the Vela neutron star afte r the 1985 glitch. 
The error bars are comparable to the width of the line; see £14.51 



(|52"j) must be positive definite. Choosing K — 0.001 and p n — 0.1, we deduce 
I/J = 6.898 and (3E 1 / 2 = 2.75 x 10~ 6 from ([!)]) -([M]). Using Figured we find 
B = 4.20 x 10" 6 and B' = -1.05. For these values, we calculate J = 0.120, 
and hence, from ([3"0")) . we obtain E = 1.05 x 10~ 5 . The extracted value of E is 
roughly consistent with the theoretical expectations outlined in the previous 
paragraph, however B and B' are not. The agreement is unsurprising in view 
of the simplifying assumptions made, e.g., cylindrical geometry and initial co- 
rotation between the fluid components. Glitch recovery in neutron stars has 
been explored more fully in a recent paper, where an analytic solution in a 
sphere is applied [30]. However, because neutron stars are strongly stratified, 
the extracted parameters must be treated as body averaged quantities, an 
issue which does not affect terrestrial experiments. 
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5 Conclusions 

We present a recipe for interpreting non-invasively the transport processes at 
work in supcrfiuid spin-up experiments like those pioneered by the Tsakadze 
group |15j , whose observed features have resisted theoretical explanation until 
recently. The interpretative recipe is sensitive to certain physics, e.g., Ekman 
pumping and supcrfiuid turbulence, that does not occur normally in second- 
sound or torsional-oscillator experiments, which rely on small amplitude vi- 
brational modes. By analyzing the response of a supcrfluid-fillcd cylinder to 
an impulsive initial acceleration, the recipe also yields unique values for the 
transport coefficients B, B' and rj in general, two relations between the three 
coefficients in the strong coupling regime, and can be extended to vessels of 
arbitrary shape |311I32| . Indeed, in the strong coupling regime, two identical 
spin-up experiments in differently shaped containers can be enough to lift the 
degeneracy between B, B' and E inherent in a single experiment (see section 
£13.31 and £14. The method cannot compete with traditional measurements 
of B, B' and E in He II but may help supplement other techniques in 3 He 
and 3 He- 4 He mixtures, where accuracies at present are lower. The method 
is a powerful tool in interpreting the spherical experiments of Tsakadze and 
Tsakadze [109U110| fT5 ] , and glitch recoveries in neutron stars [30] , where only 
remote observations are possible. 

CAVE acknowledges the financial support of an Australian Postgraduate 
Award and the Albert Shimmins write-up award. CAVE also thanks the 
anonymous referees for their constructive suggestions on the manuscript. 



Appendix 



For a container of arbitrary axisymmetric shape, the fluid is contained within the 
volume —h(r) < z < h(r), where (r, z) refe r to cylindrical coordinates, scaled by the 
length-scale of the container L. Equation (|13|l holds, with the modified definitions 



K = I}/Ic, 



a (t) 

9 (t) 

P 
I(r) 



An 



dr r 3 J(r) [ e <"+M r - e u -^ T ] 



u+(r) -u-(r) 
47r f R drr 3 {to-(r) [ e "+ (r)r - l] - u+(r) [e"" (r)T - l] } 



'/ Jo 



U+(r) - u)-{r) 



h(r) 



± i 



h(r) 



h(r) 



1/2 



2BE~ 1/2 
2 - B' 



A_(r) [Xl(r) + Xl(r)] 



1/2 



[A^(r)-A 2 _(r)] 2 +4A^(r)A 2 _(r) 



(35) 
(36) 

(37) 
(38) 
(39) 

(40) 
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J(r) = 



A±(r) 



Pn H 2 (r) 
2pA_(r) [l+# 4 (r)] 

x{[A 2 (r) + A 2 + (r)] [l - iJ 4 (r)] + 4A 2 _ (r)i/ 4 (r)} 



£T(r) 



(B'-2) 2 +B 2 



ff(r) = 



If = 



_(p n B'-2p) 2 + (p n B) 2 

2p s g [l + g»] 
H 2 (r) [{p n B' - 2pf + (p n B) 2 

1/4 



1/2 



1/2 



1 + 



pL 5 



= 4tt 



drr 3 h(r 



(41) 



(42) 

(43) 
(44) 



where If is the scaled moment of inertia of the contained fluid (as if it were rotating 
as a rigid body) and R is the cylindrical radius of the ves sel, s cale d to i. In a 
cylinder, one has h(r) = 1, L = h, R = r/h, and we recover (fT4^) - (f2T|) . 
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